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Abstract. We study the Sinai model for the diffusion of a particle in a one dimension random 
potential in presence of a small concentration p of perfect absorbers using the asymptotically exact 
f5 real space renormalization method. We compute the survival probability, the averaged diffusion front 

^ and return probability, the two particle meeting probability, the distribution of total distance traveled 

before absorption and the averaged Green's function of the associated Schrodinger operator. Our 
I work confirms some recent results of Texier and Hagendorf obtained by Dyson-Schmidt methods, and 

extends them to other observables and in presence of a drift. In particular the power law density of 
states is found to hold in all cases. Irrespective of the drift, the asymptotic rescaled diffusion front 
^ of surviving particles is found to be a symmetric step distribution, uniform for < ^^(t), where 

^ 5(t) is a new, survival length scale (^(t) = TXnt/yTp in the absence of drift). Survival outside this 

I sharp region is found to decay with a larger exponent, continuously varying with the rescaled distance 

x/^{t). A simple physical picture based on a saddle point is given, and universality is discussed. 
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1. Introduction and model 

The Sinai model U of a particle undergoing thermally activated diffusion in a one dimensional 
random energy landscape U {x) is the simplest tractable model of glassy dynamics and has generated 
considerable attention in physics, in probability theory and even in biophysics and in finance [2|3JS1E]- 
In this model U{x) has the statistics of an unbiased random walk, i.e. {U{x) — U{x'))'^ ~ 2a\x — x'\ 
hence the barriers grow as \/\x\ with spatial separation and as a consequence diffusion is ultra-slow as 
X ^ hi' t. Interestingly, the associated Fokker-Planck operator maps to the Schrodinger operator of 
a ID supersymmetric quantum mechanics describing a class of particle-hole symmetric random hopping 
Hamiltonian. These have been studied using replica [2], Dyson-Schmidt [BJ HJ [7] and supersymmetry 
methods [8J. They exhibit the Dyson singularity N{E) ~ l/(lni<^)^ of the (integrated) DOS and the 
exact (quasi-delocalized) eigenstate at E — Q and are ID analogs of the heavily studied 2D chiral class 
which in two dimension exhibit delocalization phenomena. The Sinai model enjoys many remarkable 
properties, among them anomalous drift x ^ t^^ in presence of a bias in U{x), followed by a transition 
to non-zero velocity [3]. Another remarkable property [TU], at zero bias, is that, up to rare events, the 
thermal packet is concentrated in a finite region around the bottom of the potential well available at 
time t. This property has allowed to apply a powerful real space renormalization method (RSRG) to 
this problem and derive a number of exact results [11] . The idea is to eliminate iteratively all barriers 
smaller than F = Tint and leads to results asymptotically exact in the large F limit. This method was 
introduced in the closely related context of random quantum spin chains where it has proved extremely 

useful [niin]. 

Recently, Texier and Hagendorf (TH) have studied the Sinai model in presence of dilute perfect 
absorbers and found |14j an interesting power law decay of the average probability of return to the 
origin, hence a power law behaviour of the density of states. This result was obtained using Dyson- 
Schmidt and related methods [3 [T5]. The aim of this paper is to study this problem using the RSRG 
method. Our results confirm some of the ones of TH and in addition we compute a number of other 
observables. 

Let us consider the Sinai model, i.e a random walker performing Arrhenius diffusion at temperature 
T in the landscape U{x) = V{x) — fx, i.e we allow for a bias towards the right / > and define 
5 = f /{^(j) > 0. We set {V{x) — V{x'))'^ ~ 2\x — x'\, hence choose units of length such that 
(7 = 1. We do not here redefine the various realizations (discrete or continuum time and space) 
of the diffusion model, and we refer for that to Section II of Ref. [IT], our parameters being identical. 
If we now add a Poissonian absorption potential W(x) with infinite strength, i.e a Poisson distribution 
of perfectly absorbing obstacles, the system is in effect cut into independent segments, i.e intervals 
between scatterers. Let us consider a model with impurity density p and distribution of the length 
of segments between impurities is p{L) = pe~P^ . We are studying the dilute limit where the average 
distance Lp = 1/ p is much larger than all other crossover or microscopic lengths. We assume no global 
constraint on the energy landscape U{x), i.e. it performs a free random walk (i.e. in the universality 
class of the free Brownian motion at large scale, upon rescaling). In that case the landscape in each 
segment is also a free random walk (which assumes, however, the same value at the common boundaries 
of each pair of consecutive segments) . 

The outline of the paper is as follows. In Section [2] we recall the RSRG method for the standard 
Sinai model in system of finite size L with absorbing boundaries. In Section [3] we average over the 
distances between absorbers and obtain the survival probability and its power law decay in time. In 
Section [4] we compute the full averaged diffusion front and analyze its form at large time. The survival 
length scale ^ {t) is introduced. In Section [s] we obtain the average return probability. In Section [6] 
we focus on the distribution of the relaxation rates, and obtain a simple physical picture of the power 
law in time relaxation in this model. We also analyze some rare events, and whether they could lead 
to slower decay. In Section [7] we consider the diffusion of two particle and their meeting time. In 
Section [8] we compute the distribution of distance traveled before absorption. In Sectionj9]we obtain 
the Schrodinger Green's function, and a discussion and conclusion is given in Section |10[ Appendice 
A contains more on two particule survival and Appendix B on the diffusion front in the absence of a 
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Figure 1. Schematic representation of the finite size RSRG measure[l]for the renormalized landscape 
{x axis is length, the £i are the bond length, and the vertical y axis - not drawn - is the disorder 
potential energy). It is a sum of terms with fc = 1, 2, .. valleys (i.e. sets of two consecutive bonds) 
not yet decimated at scale F = Tint, i.e. of barriers larger than F (in vertical y direction). The 
edge bonds (first and last) cannot be decimated, only the bulk bonds (all the others) are. In the 
absorber model, there is an absorber at a; = and x = L. A particle which started at t = on the 
left edge bond (resp. right) has been absorbed at time t by the left absorber at a; = (resp. the right 
absorber at x = L), asymptotically with probability one. A particle which started at xq on bond 2n, 
1 > n > A: has survived and moved to the right at x, the bottom of the starting bond (asymptotically 
with probability one), while a particle starting on bond 2n + 1 - not shown - has moved to the left 
to the same point x (asymptotically with probability one). 



random potential. 

2. Finite size measure 

Let us recall first the finite-size RSRG for the Sinai model. More details can be found in [TT] and in 
e.g. [inill^ for quantum models. The measure for the renormalized landscape at scale F for a system 
of fixed size L can be expressed as a sum: 

NrM - lrE^{h)E+{£2)S{L - {£, + £2)) (1) 

00 k 2fc+2 

+ J2lrE^{£i){l[P+{l2,)P-{e2,+i))E+{£2k+2)S{L~ ^ £,) (2) 

k=l j = l 1=1 

i.e. a sum of measures for the events where there remain 2fc + 2 bonds in the system, i.e. k valleys 
of bulk bonds and two edge bonds, with fc = 0, 1, ... We have singled out the term k = which plays 
a special role. The product form reflects the Markovian nature of the landscape, which is preserved 
by decimation, hence is either an exact consequence of the choice of a Markovian initial landscape, or 
a consequence of the convergence to the fixed point landscape. The only constraint is the fixed total 
length, implemented by the delta functions. Bonds in the bulk have a length distribution P^{£) while 
edge bonds have length distribution E^{£). The subscript + denote descending bonds along the bias, 
while — are barriers which oppose the bias. As seen in Fig. 1 the left edge bond is ascending, and 
corresponds to the " absorption zone" of the left boundary, while the descending right edge corresponds 
to the " absorption zone" of the right boundary. This geometry correspond the finite size measure with 
so-called AA boundary conditions, i.e. absorbing on both ends as relevant for the present problem, also 
shown in Fig. [2] The factor Tp which is nothing but the average bond length ensures the normalization 



Sinai model in presence of dilute absorbers 



4 



of the total probability to unity. Let us call this probability. The measure is easier to write in 
Laplace transform with respect to L, i.e. multiplying ([T|) by e~^^ and integrating dL one gets the 
normalization condition dLZ^e^P^ = 1/p, i.e. = 1, in the form: 



l-P+{p)Pf{p) p 



It -rvi^-rv/- i (3) 



The fixed point form for the bond length probability, in Laplace Pr{p) — Jq e P^Pr(^), takes the form 



y/p + 6^ cosh(rVP + <52) T Ssmh{T^yp + S^) 



j^±, ^ ^""^ T /Sinh r(5 ,2 

smh((5r)(^p + (52 coth(r^p + ^2) ^ ^) 



which is easily checked to satisfy (|3|. 

Note the simpler form in the absence of bias: 

Pr(p) = ^^ , Erip)^^!^ , Ir^T^ (6) 

cosh(r^) 

and the expressions after inverse Laplace transform (denoting for simplicity the function and its LT 
by the same symbol): 

-|-QQ -|-00 



-oo n— — oo 



We are now equipped to study the problem with absorbers. The idea is to use the efi^ective 
dynamics: a particle starting at time t = at xq is at time t at the bottom of the bond (or the 
valley) containing xq in the renormalized landscape at F = Tint It is there with a probability 
which tends to one as F ^ oo (corrections due e.g. to barrier or well degeneracies, or anomalously 
rare thermal and disorder configurations can also be estimated). In the absence of absorbers this 
leads to a diffusion length scale given by the average bond length, L{t) = Ir = T'^ = T^ln^t in the 
symmetric case. Let us stress that while the RSRG is valid for any 6 to determine landscape extrema, 
its (simplest) application to the Arrhenius diffusion requires that the bias S is small. Hence the results 
presented here will be valid to lowest leading order in 6. Extensions to higher orders in S are possible, 
but more difficult as the effective dynamics must then be modified to also take into account multiple 
well contributions order by order in S, see Ref. [TH]. Note finally that we denote disorder averages by 
— but use rather the notation (..) ^ for specific averages over the finite size block measure ([ij, although 
both have the same origin, i.e. the random potential U (x) (we do not need here a notation for thermal 
averages) . 

3. Survival probability 

Let us first compute the survival probability S(t) of a walker up to time t. We assume that the walker 
starts at xq with uniform measure on space. Hence the probability density that it starts on a segment 
of length L is Lp{L)/ J dllp{l) — pLp{L). Given that, it has a uniform measure on the segment and 
its probability Si,{T) to be alive at time F = Tint is equal to the probability that a;o is not in the 
absorbing zone, see Figs. [T]and[2] equivalently that is on a bulk bond: 

(^1 + hk+2)L ^ {h + ^3 + ■■ + hkjL j-g-j 
L L 



It is easy to compute the Laplace transform, using the normahzation identity ([s]): 

/dLL^.(F)e--^l + |;^ + ^ (9) 
J p2 pE+{p) pE^{p) 
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Figure 2. Absorbing boundary in the RSRG. Top: the initial landscape U(x) for x > 0: Uq is 
artificially set to a very large negative value. Bottom: the renormalized landscape after the iterative 
decimation of the smallest barriers up to barrier size F. The fist edge bond - which is never decimated 
- is an absorbing zone: any particle starting there at time t = has moved to the left edge while 
encountering only barriers smaller than F, hence it is absorbed with probability asymptotically equal 
to one at time Tint = F. Particles starting on the next (i.e. bulk) bond have moved to the right 
and have survived. Another such boundary is present at x = L (not shown here). 



Setting p ^ p this immediately gives the survival probability in the original model, averaged over all 
segments. Let us first analyze the symmetric case: 

S{t) = p- I dLLSdT)e-^^ = . = . l':' , (10) 

^ ' ^ J ^ ' sinh(2r^) sinh(2yprint) ^ ' 

= l-lpT^lnH + 0{{Tlnt)^) , T\nt<^l/^ (11) 

«4V^ (Tlni) t-2Typ ^ rint>l/y^ (12) 

in respectively the small and large time limit. S{t) correctly decreases from 1 to zero as a power law at 
large time, plus logarithmic corrections. The crossover time occurs when the typical difl^usion length 
in the absence of absorbers L{t) — In' t becomes of the order of the typical separation between 
absorbers Lp = 1/p. The validity of these results is small p and large times, i.e. both lengths L{t) and 
Lp sufficiently large so that the renormalized landscape has reached its RSRG fixed point form. This 
implies that diffusion has crossed over from pure to ultra-slow (so short times here are always meant 
as large compared to the crossover time from pure diffusion to ultra-slow [21 ). 
In presence of an applied force, i.e. a biased landscape one finds: 

Sit) ^ p- I dLLSUr)e-^- = /' + pTlntvVT^coth(Tln^7;;T^) 
7 2d^ + p+pcosH2Tliity/^TP) 

a generalization of the previous result. The general validity of this result is now small p ^ S^. Let 
us recall that in the absence of absorbers there is a crossover length scale Lg = 1/(5' below which 
the diffusion remains ultraslow and almost unaffected by the bias. Hence it is not surprising that 
the condition for formula (13 1 to crossover to the result for the symmetric case (10 1 is ^ p, i.e. 
Lp < Lg. 

We can now define a small time regime Tint < 1/y/p+S'^ where: 

Sit) = 1 - ^pT' In' t - ^/o(8<5' - 7p)T^ In* t + 0((rint)S) (14) 
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and a large time regime Tlnt^ Xj ^ p + P- where the formula ( 13 1 can be approximated by 



^(t).4m±Z:^^ilv^ , Tlnt»l/y^ (15) 
In the general case p ^ 5^ this yields the large time decay of the survival probability: 

X2 

S{t)K^{—+T\nty/p+5-^) t-^T^p+s^ ^ TlnO l/v/p + <52 , p 5^ (16) 
P 



a simple generalization of the symmetric result (12 1. In the case where p ^ (5 , i.e. Lp '3> Lg 



there is however another crossover time in the large time regime. As clear from (151, it is defined by 
1 <C \/p+ 6'^Tlnt ~ 41n((5^/p), i.e. time around t ~ (J^/p)''/''^ in which case: 

S{t)^ 1/(1 + ^t'^^) (17) 

crosses over from S(t) = 1 to a small S{t). This additional crossover time is necessary since one 



must recover a unit survival probability when p — s- 0, and that is not obvious on (16). Hence when 
p <C i5^, i.e. Lp ^ Lg, the large time regime splits in two subregimes depending on whether (i) 
1 < a/p+ S^Tlnt < 41n(J^/p) in which case S{t) w 1 and (ii) ^/p+ S^Tlnt > 41n(5^/p) in which 
case S{t) is again given by ( [T6| . 

The decay is thus faster in presence of a bias. This is intuitively clear since in the absence of 
absorbers the motion of the particle is faster at large time and becomes x{t) ~ t^^'^ /(AS^): hence 
the particle should die sooner. This effect however does not make the decay of the average survival 
probability exponential (or stretched exponential as in the absence of the random potential) , here it 
remains a power law in time. This is because of the rare regions associated to a barrier against the bias 
which contain particles trapped for a long time. As found in ^llj in presence of a small bias 6 1 and 
without absorbers, the diffusion front at large time takes the form P{x, t\0, 0) « 9{x) exp{—x/x{t)) /x{t) 
hence the probability that a particle has remained in a region of size ~ 1/p near its starting point is 



l/(pa;(t)) ~ AS'^t "^^"^ / p, a decay which exhibits some similarity with (16 1 taken for p <C 5^. In the 



other limit p':S> S the effect of the bias is weak and one recovers the result (|12| for the unbiased case. 
4. Diffusion front 

Let us now compute the fully averaged diffusion front P{x, t\xo, 0), i.e. the probability that the particle 
has survived up to time t and has moved from xq to x in presence of the absorbers. 

One first defines PL{x,t\xo,0) in a fixed segment of size L for the problem without absorbers, 
where here the average is only on U{x) and one studies its Laplace transform P(p,pQ,q) = 
dL dxo dxe~P^~P°^°~'^^ Pl{x, t\xQ, 0). Since in the effective dynamics the particle is at time 
r = Tint at the bottom x of the bond containing xq, for each term containing fc > 1 bulk valleys in 
the measure ([ij one first defines: 

2n 2n-l 

Pik+2 - {H^ - Y.^^^^0< x))L,k (18) 

1=1 i=l 
2n 2n+l 

P2kt2 - (^(^- - E ^^)^(^ < ^0 < E (19) 
i=l 1=1 

corresponding to the initial condition xo being either on the 2n-th bond or the 2n + 1-th, while x is 
at the bottom of the 2n-th bond, with n = 1,2, ...fc (see Fig. [iJ. Note that the P^" corresponds to 
the particle moving to the right, i.e. y = x — xq > while the p2"+i correspond to motion to the 
left y <0. Hence one separates Pl{x, ^I^Oi 0) = ^ti^^ A^o, 0) + P£{x, t\xo, 0) into the contribution of 
right motion (P"*") and left (P~). In Laplace one has: 

-F'2fe+2(P'?'0: 9) ~ ^^p+po+qiPp+Po+qPp+Po+q'^ (Pp+q ^ Pp+Po+Q^ ) ^q 
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Using the sum rule (|3| and summing over n — 1 and k — n one finds: 

PHp,P,,q) = I , \ ^ iP:+q-Pp\po+q)P,- (20) 

f rPoq[P + PO+q) Ep+q+p„ Eq 



J (Pq Ppo+q^Pp ^ 

irPoq{p + PQ + q) E+_^g^p^Eq 



where P^{p,po,q) are the Laplace transforms of (x,t\xoO), corresponding to motion to the right 
y > and to the left y < 0, respectively. 

We can now obtain the fully averaged diffusion front in our original model in the form: 



P(xo + y, t\xQO) = e{y)P+{xo + y, t\xoQ) + e{-y)P-{xo + y, t\xoQ) (22) 

as the sum of contributions of the P^" and p^n+i probabilities, respectively. The first one is obtained 
by setting = —p, i.e. one has: 



dye~Pyp+ixo + y,t\xoO) (23) 
= dLLe-P^L-^ dxo dye-Pypl{xa + y,t\xaO) (24) 

Jo Jo J-xa 

^ '^''^^ £rP E+Ep 



while the second one is: 



dyePyp-ixo + y,t\xoO) (25) 

L-Xq 



= pM dLLe-P^L-^ I dxo dyeP^ P£ {xq + y,t\xoO) (26) 

Jo Jo J-XQ 

= P'^P~{-P,P,P) = 



r.'2u~t ^ „ ^\ — ^ Pp + p)Pp^ 



irp E+Ep 



4-.1. symmetric Brownian landscape 

Let us first analyze this result for the symmetric case: 

p cosh(r\/p) 

p sinh2(r^) 'cosh(r^) cosh(rVp + p) 



dye-Pyp+ix, + y,t\xoO) = ^ T^J^Z i ^.^.n. ^^ - .^^ ) (27) 



r^p i + rvp(2tanh(ryp)-coth(rvp)) ^ ^^^^^ 



. ^.-[1 '-^^p+0{p')) (28) 

smh(21 y'p) Ap 

and the same result for P^ with e'"^^ e^^ . The first term in the small p expansion is half of 
the survival probability obtained in (10), as expected. The second gives the average displacement 
conditioned to having survived (i.e. averaged over the walkers which have survived) is: 

_ 1 + rVp(2tanh(r^) - coth(r^)) 

\y\ - 4^ i^yj 

with r = Tint, with the following asymptotics respectively and small and large time: 
5 29 

= -T'lnH-—pTHiiH + OilnH) , Tlnt^l/^ (30) 
= ^Tlnt+^-l^t-'^^P+.. , Tlnt»l/^p (31) 
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The small time behaviour \y \ — In^ t is the one obtained in the Sinai model in absence of absorbers 
[TT] . Note that the distance traveled by the surviving walkers only grows at large time as ~ ^(t), with: 

at) = Tint/ ^ (32) 

as these walkers are those which remain closer to their initial position (within a region whose size is 
still growing with time). This new length scale ^(t) is discussed further below. 

Performing the inverse Laplace transform from (27 1, one obtains the fully averaged diffusion front: 

21 



7rpcosh(r^) ^ (-l)"(2n + l) ^.2'^^\y 



+ ^' ^l-'^) ^ sinh-(r^) h Pr^ + ^in+^^ - + 2)^T^) 

with r = Tint. It is independent of the starting point and normalized to the survival probability 
dyPixQ + j/,i|xoO) = S{t). In the limit Tint ^ ^/y/P where the probability that the walker has 
met an absorber is still small, one recovers the Kesten diffusion front derived in [11 using the 
RSRG (the same as the limit p — Q). 
Using the Poisson formula: 

oo oo Z'+OO 

E E / dxf{x)e'^-^^ (34) 

n— — oo m— — oo ^ 

This can be rewritten in the dual form as: 



sinh(rVp) ^ TT 



\y\ 



At large time T ^ Tint ^ l/x/P' ^^^^ dominated by the term to = 0, which yields: 
P{xo + y, t\xo, 0) = S{t)^e^^P j^^ ^^""^"^ 



Now one can perform the change of variable y — Ty/ ^fp and z = Tz/ ^fp in (35 1 and taking the large 
r limit one easily sees that either the saddle point point at 5 = 1/2 is in the integration domain, and 
the result is constant, or it is not, in which case it gives a subdominant contribution. Hence in the 
large time limit the diffusion front for the rescaled variable y = y/£,(t) (where ^t) is the new length 
scale defined above in ( [32^ ) takes a remarkably simple form: it becomes a uniform distribution in the 
interval \y\ E [0, 1/2] and zero outside: 

Vp \y\^. 



P(xo + 2/,t|a:o,0)«5(i) ^ 0(--^) , rint»l/V^ (35) 

where 0(x) = 1 for x > and 0(x) — for a; < is the Heaviside step function. The first moment 
of this distribution reproduces (31 1 above. Note that each integral in (35 1 has a saddle point at 
z = {2m + l)r/(2Y^), resulting also in a step function probability for \y\ up to the saddle point, but 
the higher values to > 1 contribute only as a subdominant relative factor t~2r?xTyp ^ S{t)"^ to the 
total probability, as compared to the main contribution. 

Equivalent manipulations can be performed on the Laplace transform. The limit for large F of 
the expression (27 1 reads: 

' e-Pyp{x^ + yAxo,Q) = S{t)^{l - e'^iVP-^p)) 

pT 

where we recall that LTy^p dzf{z) = ^{f{p = 0) — f{p))- In this expression further approximating 
y^p + p = y/p + 2^ + -. one recovers the Laplace transform of (psl) on the positive side. 

To summarize, inside the sharp region (survival region) — ^ < y — y/^t) < ^ the diffusion 
front is uniform and decays uniformly as S{t). Outside this region, there is a faster decay, which, 
in the variable y can be described as power law in time with a continuously varying exponent, i.e. 
P{xo + y, t\xQ, 0) with a decay exponent: 

a^T^p{y+^ + l)>2T^p (36) 
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for fixed y > 1/2. 

It is interesting to note that for free diffusion, i.e. no random potential U{x) ~ 0, an analogous 
phenomena occurs, i.e. there is also a surviving zone — Co(i) < 2/ < S,n{t) with ^o(^) ~ t^^"^. It is sharp 
but not as sharp as what is found here for the Sinai diffusion. The details are recalled in Appendix B, 
where the analogies are emphasized. 



4-. 2. biased landscape 

The formula for the biased landscape obtained by inserting Q into our general result ( 23|25 1 are more 
cumbersome, so we will not display them here, but concentrate on the large time limit. 

Let us first display the normalizations, i.e. the survival joint probabilities that the walker has 
survived and moved to the right P+ (resp. and moved to the left, P~): 

S^T±S + TpT 6Ty/S^ + pcoth (Vy^WTp) 



S^{t)= / dyP±{xo + y,t\xoO) 



with r 



V'PTpsinh (2Ty/S^ + p'j ±6{1~ cosh (2rV^ 
Tint. This yields the asymptotic behaviours respectively at small and large time: 



17Sp 



T^ln^t + 0(ln*t) 



= 2(TlniV7+^ ± 



y'S^ + pTS' 

These two probabilities correctly sum up to the total probablity S{t) 



-P) 

(37) 
(38) 



S+{t) + S-{t) that the 



walker has survived (|13l 
Tint 



5^ and the crossover then occurs at 



The validity is again small p 
\ I \l~p~-\~P' ■ In the case p <^ 5^ the formula (38 1 is valid only for much larger times, as 
discussed in section |3l 

Remarkably, the large time behaviour ( 38 1 shows that the dominant term is the same, i.e. surviving 
particles have the same probability to have moved to the right than to the left. Let us recall that for 
(0 = one finds: 



strait) — Pright — 



1 



coth((5rint) 



6Tlnt 



(39) 



2sinh^(JTlnt) 

S-^^{t) = Pie ft = 1 - Pr^ght (40) 

Hence Pieft ~ {25Tlnt — l)t^^^^ and at large time almost all particules go to the right. One 
notices that for < p <C (5^ one has S-{t) « Pie ft hence in that regime the particles which go 
to the left are assured to survive, while the one which go to the right are assured to die. More 
generally, the ratio 5+(i)/S'^Q(i) sa S^{t) ~ i-'^T\/p+s^ large time, while for left going particles 

it is s-{t)/s;^,{t) « 

i.e it has a slower decay. Thus survival is always better among 
the particules which go to the left, but these are a small minority of them as the drift pushes most of 
them to the right. At the end the number of survivors is the same in the two vastly unequal right and 
left-moving populations. 

Examination of the large time limit of the Laplace transform of yields a result very similar 
to the symmetric case, i.e. formula (36 1 replacing p by p + (5^ everywhere. There is a surviving scale 
which reads now ^ = {Tint) / \/ p + 6'^ . Expanding further in the scaling region of fixed pT, i.e. fixed 
y ~ y/ii inverting the Laplace transform one finds again: 

'2 



P{xo + y,t\xoO)~ Sit) 



) , Tint > l/v/p+^2 



(41) 



Tint '2 Tint 

i.e a fully symmetric distribution for the surviving particles, with the same step function shape. Again 
this is valid when p ^ S^. A more complicated crossover occurs in the large time regime when p <C (5 as 
the diffusion front first takes the forward exponential form ^ 9(x)exp(x/x{t))/x{t) as in the absence 
of absorbers and finally crosses over to (41 1 at a much larger time which diverges as p 0. We will 
not study this crossover in details here. 
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5. Average return probability 

A slightly simpler quantity is the Laplace transform of the return to the origin P(po,q) = 
dL dxoe~P''^°~'''^ Pl{xo, t\xo, 0). It is easily seen that one just needs to sum: 

Pik+2iPo, q) = ^"r£;p-+,(P++,PpT,+,)"-^Pp++,P-(P+P-)'=-"i?+ (42) 

leading to: 

P{Po,q)= . ,\ (43) 

The fully averaged return probability in the original model is given by summing over all initial 
positions with uniform weight and then over all segment sizes giving: 

P{xQ,t\xo,Q) = / dLLe-P^L^^ / dxQPL{xQ,t\xo,Q) = - , ^_ - p (44) 



irEp Ep 



This gives in the general case: 



P(xo, i|xo, 0) = / , = , (45) 

sinh2(ry;^+J2) sinh'(vVT*2rint) 



and the large time asymptotics: 



P(:Eo,i|a;o,0) -4(p+52) ^-2tV^ ^ rint< + (45) 

In the case of a symmetric landscape (zero bias) this is in agreement with the results of TH [14] (to 
compare one should set T = 1 here and (7=1/2 there) in the limit of small p (to leading order in p) . 
Here we obtain the averaged return probability in the biased case as well, which exhibits a very similar 
expression. In particular, contrarily to other quantities discussed above it exhibits no complicated 



crossover for p ^ 5^ . Indeed one can directly take the limit p ^ in the large time asymptotics ( 46 1 
and recover the probability of return to the origin of the biased Sinai random walk which is also a 
power law [2j[9l[IT]. 



6. Relaxation rates and physical picture of the power law decay 

The aim of this Section is to explore in more details the distribution of relaxation times and obtain a 
simple physical picture for the power law relaxation. 



6.1. lowest energy level 

The (left and right) eigenfunctions of the diffusion (i.e. Fokker-Planck) operator (t>L.R{x) are related 
to the eigenfunctions of the Schrodinger operator ■0(a;) as (t)R{x) = e~^^^^/'^'^'ip{x) and (t)L{x) = 
eC^(a:)/2T^^^^^ This relation is compatible with the vanishing of both eigenfunctions at the edges 
of the segment [0,L], as required here for perfect absorbers (Dirichlet boundary conditions). The 
eigenenergies Ei, i.e. the relaxation rates, are thus the same for both operators. 

The RSRG method allows to find all lowest eigenenergies as = —TlnEi where the are the 
scale of the successive decimations [12] • By convention we number < Eq < Ei,.. the energies in 
increasing order, hence Fq > Fi > ... the corresponding decimation scales. Other methods to compute 
the statistics of the eigenenergies are presented in Ref. [20, and in Ref [7], the latter being discussed 
at the end. 

We now recall the calculation of the probability distribution of the lowest energy level Eq, i.e. 
the smallest decay rate in a given sample of size L. It is associated to the decimation of the last bulk 
valley and occurs in the block fc = 1 of Fig. [T] when F reaches the size Fq of the smallest barrier of 
the two bonds, bringing the block A: = 1 to fc = 0. Hence the probability Pro6i(Fo < F) that the last 
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decimation Tq = —TIuEq in the segment [0,L] has akeady occured at scale T, i.e. that Fq < T, is the 
total weight of the fc = term in the measure ([T|) : 



f°° - 1 

/ dLe-f'^ProbLiTo < T) = lrE^{p)E+{p) = -(1 - P^{p)P+{p)) 
Jo P 



(47) 



using the normalization From this, the probability density PL{To)dTa for Fg = — Tlni^o is 

obtained as PLi^o) — drProbL{To < F)|r=roi s.nd one finds: 

r dLe-P^PUTo) = dr„ {- — — ^-^^=— -) (48) 

Jo {p+S^)coth\ro^yp + 6^)-S^' 

Note that this formula was already obtained in [IT] (Section VII-B) for the equilibration time, i.e 
the (inverse of the) energy of the first excited state, in a finite size sample with reflecting boundary 
conditions. However, in the limit studied here of a large sample and large times this turns out to be 
the same distribution (in a rescaled sense) as the lowest decay rate for the present problem, i.e. with 
absorbing boundaries [23] 

In the absence of bias one finds: 

/ dLe-^^Pi(Fo) = -aro( , ^ = -^d,{ ) (49) 

Jo pcosh [Toy/p) lo cosh (Fq^) 

One can check that the relation ^^^^1^]^^) ^ -^r(p) + 2p9p£'r(p) holds in that case, and one finds, 
upon inverse Laplace transforming and using ([T]): 

PLiTo) = ^{LEr„{L) - 2LdL{LEr,{L))) (50) 
i 



2L ^ , 1^9 i ^ -'^^("+5)^7^ . ^ 



Fi^ ^ ^ ^ 2' 

^ n— — oo ^ 



^(-l)™+im2e-™ T (52) 



^ m— 1 

using Poisson formula. 

Let us now recall that in a given segment [0, L] the return probability is a trace over all eigenstates 
/g dxQPL{xo,t\xQQ) — X)n^~^"*- Going back to the model with absorbers one averages over all 
segment sizes and starting points and one finds the relation: 

P(xo,i|xo,0) = p2 r dLLe-P'^ r ^PL{xo,t\xoO) = p' H dLe-""^ V (53) 
Jo Jo ^ Jo „ 

This quantity was computed in Section [5] To get physical insight into what controls the large time 
power law decay we will keep only the ground state in the sum, hence we get only a lower bound on 
the total decay. However we want to estimate the thermal average more precisely, to test the accuracy 
of the RSRG (i.e. to make sure we have not overlooked a process leading to a slower decay). From 
(53 1 we thus obtain in the absence of bias 6 = 0: 

P{xo,t\xo,0)>p^ / dLe-P^ dToe-'^'""" Pl{Vo) = P dT^e-'-'""'" d^.t^nh^To^p) 
Jo Jo Jo 

= Apf^ dy^l_|_e-*«^'<^^^' « 4pF(l + 2ryp)t-2^VP (54) 

setting y = e~^^^^ and the last approximation is at large time. It is interesting to note that this is 
exactly the same prefactor for the large time decay as obtained by TH [14j. The effective dynamics in 
the RSRG instead replaces e~'^ function which then leads to: 

/•oo 

-P(a^oikoO) > p / dTodv, Unh^V^^p) = — « 4^^^^^ (55) 

J Tint cosh (J Int^) 
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and the last approximation is at large time. As discussed above the RSRG is only expected to be 



accurate for large Lp i.e. small p, in which case the two asymptotic behaviours (54 1 and (551 cannot 



be distinguished. Taken together the results (54 1 and (55) compared to (46 1 (in the absence of bias) 
show that (i) the asymptotic time decay can be correctly obtained by retaining only the longest 
relaxation time in each segment between absorbers and then averaging over the segments (ii) that the 
particles which jump anomalously late are sufficiently rare not to affect the leading asymptotics at 
small p [ig 

6.2. saddle point picture, universality and rare events 

We can now clarify the physical origin of the power law decay. First we can concentrate on the 
longest relaxation time (i.e. smallest energy) in each segment between absorbers to — e^°^'^ . We must 
explain how a power law distribution of relaxation times of the form dtg /^q^^"^^ is generated by the 



combination of absorbers and disorder. We note from (50) that the distribution of barriers Fq decays 
as ~ exp(— Fq/L) for large barriers, i.e. the distribution of relaxation time at fixed L has a log-normal 
tail. Averaging over the various sizes L of the segments between absorbers indeed produces the power 
law spectrum: 

" dLe-^^e'^o/i ^ e-2Vpr„ ^ ^^^t^p (56) 

via a saddle point a,t L = ^(tg) = Fp/y^ = (Tlnio)/^/^ which is the optimal size of segments 
contributing to relaxation times Iq. Note that this is also the survival length scale associated to time 
to already discussed above, and four times the one obtained in ( [31] ), which is quite consistent. 

These considerations also allow to discuss the universality of the present results, beyond the 
Brownian energy landscape. In problems of large time decay in random media involving some saddle 
point argument, one is always at risk of neglecting a slower process due to very rare events. By using 
here the RSRG fixed point based on the Brownian landscape we assume somehow that the saddle 
point indeed occurs in the scaling region Fq ~ \/L where the probability tail e^^"^^ is valid (i.e. large 



but still typical barriers). In fact, from the saddle point in (56) we need more, we need that this 
tail holds where the saddle point occurs, i.e. at Fq ~ ^/pL (barriers larger than typical). Hence it 
is safe to expect that it holds only for small p. To test the validity for larger p, one may consider 
e.g. the most extreme events where the landscape U{x) between two scatterers has an unusually low 
minimum with a barrier to exit of order ^ —L/2. One can evaluate the probability of such an event 
to ~ e~^^'^^ (e.g. for a ±1 discrete random walk U{x)), and the total contribution of all such events 
as '--^ / dLexp{—{p + ln2)L — t/e^/'^'^) ^ ^-2T(p+in2) xhis is again a power law, which turns out 



to be faster than the power law ( 56 1 for any p. Hence the Brownian landscape is stable at least to 
such events for a discrete bare model. Generally these are events where barriers are smaller than 
expected from a Brownian landscape (e.g. there is a maximum possible barrier — io/2 in a sample of 
size Lq). The opposite case, where there can be bigger barriers than for the Brownian landscape, may 
be problematic. For instance, although the RSRG fixed point is also valid for power law distributions 
of barriers in the original model with fast enough decaying tails [12], it is not clear whether, in that 
case, these tails may not always lead to a slower decay. We leave this question to future investigation. 

Another extension is the case where p{L), the distribution of segment sizes, is not simply an 
exponential. Although exact calculations can be performed using the above methods, we can sketch 



the result of the saddle point argument. From (56) one sees that if p{L) is delta-peaked on a given Lq 



the large time decay will be as ~ exp(— (T~ln t)/Lo). But then, for a discrete model, it must change 
at very large time, i.e. Tint ^ Lq, to a faster g-^-o in2-t/e^o decay as the distribution of relaxation 
times hits its large time cutoff (i.e. the largest possible barrier in a sample of size Lq). If now the 
tail is p{L) ~ exp(— (pL)'') with 7 > 1 one expects a faster decay than in the exponential case, i.e. as 
~ exp(— c-y(y/print)^^/*^'''+-'^^) via a similar saddle point argument, and ~ (1 +7)7^^/*^'*'+"'^^ which 
reproduces the previous result for the Poisson case 7=1. Again one must check a posteriori that 
the Brownian landscape tail e~^o/^ holds for the given model around the saddle point L ~ Fq^^^'*''''''. 
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And one finds, as before, that for a discrete model, one can hope that the saddle point behavior holds 
only up to Tint ^ (l/p)'*'/^'^^^) when the maximal barrier is reached. For the Poisson case 7=1, 
the saddle point picture nicely holds for all large times. Finally the above stretched-log decay formula 
should also holds for 7 < 1 when the decay becomes slower than a power law. 



6.3. discussion of the energy spectrum 

Before we close this section, let us make a few comments. First, one can also compute from the RSRG 
the distribution of the fc-th eigenvalue Ek, from the weight of the fc-th block in the measure ([T]): 

r dLe-P'^ProbUTk < F) = lrE^{p)E+{p) Y,{Pr (p) Pr (p))' (57) 

= ^(l~(Pr-(p)Pp+(p))'^+i) (58) 

In the symmetric case using Tdr = 2p9p acting on the fixed point probability, this gives: 

Pl(F,) = +^(L) (59) 

i.e. L times the 2k + 2-th convolution of a bulk bond length distribution, i.e. the distribution of size 
of a renormalized valley with 2A: + 2 bonds. 

Second, let us comment on a relation between the method used by Texier in Ref. [7 and the 
RSRG. As discussed in [19] for the infinite sample, i.e. the bulk, the zeroes of the eigenstates with 
energies TlnEi < F are on hill tops of the renormalized landscape at scale F, hence the distance 
between the successive zeroes are given by the length of the valley £1 + £2. This is precisely the result 
of Ref. [7]: the P{£) there (formula (104) and (154) there) is the distribution of length of a bulk 
valley i.e. the convolution of two P{£) here (while here P{£) denote the distribution of length of a 
bond, also called A in Ref. [7]). This was obtained there from the statistics of the escape time of the 
stochastic Ricatti equation, while in the RSRG is it obtained as the fixed point of an asymptotically 
exact decimation method. Hence as far as bulk properties are concerned there seems to be perfect 



agreement between the two methods. The formula (50 1 coincides with the result of Ref . [7 (formula 
(123) and (110) there). 

These arguments must be modified to take into account the edge of the system. In the RSRG this 
is described by the edge bonds. The arguments of appear also to be correct for Dirichlet boundary 
conditions as required here. This is because the length of the system L was chosen there to be the 
interval between two zeroes of the eigenfunctions. As a result one can indeed check that the formula 



(50 1 coincides with the result of Ref . |7| (formula (123) and (110) there). More generally formula (59 1 
for the distribution of Ek coincides with formula (132) in Ref. [7]. Given this agreement between the 
two methods for calculation of the distribution of the eigenenergies with Dirichlet boundary conditions, 
it is quite natural that we find agreement for all results which can be compared. 



7. Meeting time of two particles 

The RSRG allows to compute a host of more exotic properties of the trajectories, such as aging, return 
to the origin, persistence etc.. Here we will restrict to consider the meeting properties of two particles. 

Consider now two particles distant by j/q > at time t — 0, and let us compute the probability 
M{yo, t) that at time F = Tlni (i) they are both alive AND (ii) they have already met at least once. 
This quantity is simpler to compute since it requires the two particles to have started in the same 
segment, denoted [0, L] between two adjacent absorbers (otherwise they never meet). Note, en passant, 
that the probability that two particles start in the same interval of size L is p^e~P^{l — j^)9{L < ijq), 
hence a total probability e~P'^° to start in the same interval. Even if they start in the same interval, 
they may still be absorbed before meeting. 
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As discussed in [TT] having met is equivalent to the initial positions of the particle, let us call 
them xq and xq + yo belonging to the same renormalized valley at T. The probability M{yQ,t) can 
thus be obtained from the probability ML(yo,t) that, in a sample of size L, the interval [xq,xq + yo] 
is contained in a single renormalized bulk valley (denoted 2n, 2n + I below). One has: 

oo n—k 7 2n~l 2n+l 

fc=lri=l'^0 i=l i=l 

in the notations of Section [4] given that the measure on xq is uniform in the interval. As usual one 
computes: 

M{q,p,t)^ dLLe-'^'^ e-Py"ML{yo,t) (61) 

Jo JQ 

Inserting the finite size measure ([ij), using that L = Yli='i^ performing the integration over xq and 
yo, summing over n — I and k — n and using the normalization identity, one finds, very much as in 
Section El 

M{q,p,t)^\{N{p)~N{0)-pN'm , N{p) ^ (62) 
p2 £rq^E+{q)Ej. (q) 

Upon averaging over all segment sizes one finds that p'^M{p,p, t) is the Laplace transform of M{yo, t) 
and we obtain our final result for the model with absorbers: 

Af(yo,i)-^Tp^^„^(/(p)-p/(0)) , f{p)=p E+{p)E^{p) ^ ^ 

with r = Tint and — /'(O) = S{t) is the survival probability computed in Section [s] One finds that in 
the limit p = one recovers the result of [TT] . 
In the symmetric case one finds: 

M{yo, t) = — , yo + / C dy C dzH(z) 

' sinh(2rVp) sinh2(r^)^° tanh^T^) Jo "Jo 

8p g ^4(2m + l)^yo + 2rV^(2m + 1)^(3 + 2pyo) ^ SF^ ^,(g^(„,+ 1 



p)yo 



tanh2(rVp) (4rV + ^2(2™ + 1)^)3 

with: 

1 1 +°° 2 

= 7|l7^E(-l)""^"^^-"^'^-^^ (64) 

^ n— 1 

One can check that formula ([7 1 for AI{yo,t) matches the result (126) for the meeting probability 
(denoted there 1 — Fyg(T)) in |1 Ij in the limit p = (equivalently the leading order in the limit 
p <C yo ~ r^). For small times <c yo '^-^ 1/p one finds: 

M{yo,t)^^^-y^e'-'y"/^^'e-Py" (65) 

1 ^ TT^ 

with F = Tint. One may interpret the extra factor due to the absorbers as the probability, discussed 
above, that the two particles belong to the same interval between absorbers: if they are, then the finite 
size of this interval appears to play only a sub-leading role at short times. 

Let us now study the large time limit defined as F^ ^ 1/p and F^ ^ yo with F = Tint. Then 
we can retain only the term n = 1 in (64 1 and use again the survival length scale ^(i) = T\ut/y/~p 
introduced above in (32 1 defining yo = yo/^it)- Then we obtain: 

M(yo,i)«4Fy;5(l-yo)e-2rvP + 4pF(^)i/2 / dy ^%e-i^Vp(H^) (66) 

\/P Jo Jo v^^ ' 
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For large Fy/p the last term has a saddle point at z = 1. It is subdominant for < yg < Tlni/y/p and 
in that case the large time decay is: 

M{yo,t) ^ ATliit^il ~ yo)t-^^^ ^ {1 - yo)S{t) (67) 

Hence only the amplitude is changed with respect to the (one particle) survival probability. If the two 
particle are further apart, i.e. for yo > Tlnt/y/p the leading saddle point contribution cancels exactly 



the first two terms in ( 66 1 . What remains is thus a faster decay. To study it let us start from the exact 



formula for the derivative (valid in all time regimes), easily obtained from (64 1: 

_P 

tanh2(Fyp) 

For large Fy/p one has thus: 



dy,M{yo,t)^-^£^^ I dzH{z) (68) 

yo 



/TTZ 



hence for y^ > \ we find a decay M{yo,t) ^ ^~^v^(2'o+ sjj) -v^jth an exponent which increases 
continuously with the distance. The large time behaviour is thus similar to the one found for the 
diffusion front in Section ID It can also be understood from a related calculation of the survival 
probability of two particles, performed in the Appendix A. 

8. Distance traveled before absorption 

Here we compute the distribution Q{z) of the total distance traveled, z, by a walker before it is 
absorbed. This quantity was obtained for free diffusion (i.e. in the absence of a random potential) in 
Ref. 123]. We consider a uniform density for the initial position of the walker. 

Consider a segment [0, L] between two successive impurities. To each end of the segment 
corresponds an absorbing zone of length £i, £2 respectively (see Fig. [T]), which determines if the 
walker is absorbed to the left or to the right, respectively. For a fixed L the probability distribution 
of (4, £2) is: 

gL(^i,^2)- lim erE^{ei)E+{e2)5{ii+i2-L)^e-{h)e+{i2)S{ei+£2-L) (70) 

F— >oo 

with, in Laplace and in real space, respectively: 

e±(p) = ^J=^ , e^{e)^^^S^{6V~e) + S±S (71) 

POO c\ 

0(x) = / dy—e~y" (72) 

Jx 

0(x) being the complementary error function, and e*(£) = I/^/ttI in the symmetric case 5 — Not 
surprisingly, the above probability can be interpreted in terms of extremal statistics of the (landscape) 
Brownian motion [TT1|25]. It is useful to define the full Laplace transform: 



q{p,pi,P2)^ dLe-P"^ dh d£2e-P'''-P'''''qL{£i,i2) = e-{p + pi)e+{p + p2) (73) 

JQ JQ JQ 

which is correctly normalized to q{p, 0, 0) = 1/p, i.e. J^^ 9l(^i, ^2) = 1 for each L. From this quantity 
we can obtain a number of observables for our model: 

First we obtain the probability that an interval between successive absorbers (chosen at random 
among all such intervals with equal probability, i.e. length distribution pe~P^) has absorbing zones 
(^1,^2): 

go(^i,4) = LT^,le,,p,^,MP^Pi,P2) = pe-P^'^+'^^e-{h)e+{i2) (74) 
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Second, we obtain the probability that a walker chosen at random (i.e. with uniform density 
of initial position) belongs to an interval (i.e. of length distribution p^Le^f^) with absorbing zones 
(^1,4): 

gi(4,4) - +9pJg(p,Pi,P2) = p\h+ l2)e-P'^'-+'-'^e-{h)e+{i2) (75) 

Finally let us denote z, the distance, counted with its sign, traveled up to absorption. One has 
z — —xq if Xq e [0, £i] (right absorption) and z = ii + £2 — Xq = L — xq if Xq G [^1, L] (left absorption). 
Hence one finds: 

/•oo /"OO /"OO -1 -1 

Q{z) = p2 / dLLe-P"^ / d£i / de2{-0{^z)e{\z\ < £,) + -e{z)e{\z\ < i2))qL{h,i2) 
Jo Jo Jo ^ ^ 

= 6{-z)Q-{\z\) + e{z)Q+{z) (76) 

where we used that ^ J^^ ^j^5{z + xq) = j^6{—z)9{\z\ < ii) is the probability that xq belongs to the 
first bond and that z =^ —xq. One thus has, in Laplace, Q^{pi) — ^{q{p, 0,0) — q{p,pi,0)) and 



riP2) = P^j-^{q{p, 0, 0) - q{p, 0,P2)), i.e.: 



Q±(p) = p'-{- - ,2 J 

PP {^/p + S^ ±S){^/p + p + S^ tS) 

The total probability of being absorbed to the left is = Q (0) = 5(1 ^ ^ "^^^ ) and to the right 

Q+(0) i(l + -J=), which should be seen as a scaling function of in the limit where 

both are small. Inverse Laplace transform gives: 

g±(|z|) = (fc'^l-l(l ± 1 - H6VW\)) + v/^0(v/^^)) (78) 

with 0(2:) being the complementary error function given in (72 1. In the symmetric case, (5 = 0, the 
probability of the distance traveled z, counted with its sign is: 

Q{z) = p^{^\) (79) 

normalized on the real axis z g] — cx),+oo[. Since this quantity integrates over all the history of the 
walker, the only relevant length scale is Lp = 1/p. The distance traveled by walkers absorbed up to 
time t can also be obtained by similar methods. 



9. Schrodinger Green's function 

Let us first recall the idea of the calculation in the absence of absorbers. The Green's function of the 
Schrodinger operator is (we set T = 1 in this Section) : 

Gl{x, t\xo, 0) = e^^^-^-^^-^^y^PUx, t\xo, 0) (80) 

as discussed above this extends to finite size L, consistent with the vanishing of the wavefunctions 
at the boundaries. The infinite sample case was discussed in [TT] Section VIII. Because PL{x,t\xo,0) 
behaves as exp{—U{x)) near the actual position of the particle at time t, i.e. the end of a renormalized 
bond, it is clear that the non vanishing contributions to Gl{x, i|xo, 0) will come only from (rare) bonds 
which have a second minimum degenerate with the one at the bond edge (and x, Xq near each of them, 
the order not counting, as is a symmetric function of its arguments). Hence we introduced there 
the probability Rr{£,z) that a renormahzed bond at scale T has length i and a distinct degenerate 
minimum at distance z > from the absolute minimum in the bond. Of course z must be of order 
r^, the many quasi-degenerate minima around z ~ are counted as a different (a delta- function) 
contribution (see discussion of precise definitions in Appendix E of |11]. It was shown there that: 

R^{£,z)^P^{i~z)g{z) , g{z)^^G{z/T')e-^'' (81) 
G{X) = "'^"''"^^^ - -7^2 E (1 + ^)e-'"^/^ (82) 

n=l ^ m— — 00 
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where g{z) does not depend on the direction of the bias. 

We can now compute GL{x,t\xo,0) (the average is only over U{x) in a finite segment [0, L]) by a 
simple modification of the calculation of Pl(x, t|xo, 0) in Section [i] We simply need to replace the bulk 
bond in the renormalized landscape at T = Tint containing the initial position xq by a degenerate 
bond, i.e. a R bond (if xq belongs to on an edge bond the particle is absorbed, hence we never have 



to replace those). Hence in the average (83 1 we need to write now: 

2ri 2n 

G2fc+2 = Wa;-I]^.)'^(a^-a;o-^))L,fc , G^Xl ^ i^ix -J2i^)S{xo - x - z))L,k (83) 

i=l i=l 

where x is at the bottom of the 2n-th bond and xq either on the 2n-th bond (first line) or the 2n+ 1 
(second line). Hence the average is such that in the first case the 2n-th bond is a i? bond while in the 
second case the bond 2n + 1 is a i? bond. Note that strictly speaking x and Xq can be interchanged 
(when a bond is degenerate there is an ambiguity on the definition of the two bonds forming the 
valley). We find it more convenient, and equivalent to symmetrize at the end. 



Introducing now G{p,po,q) ~ Jo° dL dxo Jq dxe G L{xt\xoO) and using that in 

Laplace i?^^ — P^gq+p, we obtain 

G2fc+2(P:P0, ?) ^ ^^^p+po+qi^p+po+q-^P+Po+q^" ^ Pp+pa+q9p+qPq ^^q ^q )^ "^g" 
^2fe+2(P'P0' 9) — ^r-Ep_,_p|j_|_^(Pp^p^+^Pj3_|_pj^_|_g)" ^Ppj^p^j^qPq gpg+q{P^Pq 

Using the sum rule (|3| and summing over n — 1 and k — n one finds: 



(84) 



ri+( \ ^ ^p+po+qdp+l^q 
G^{p,Po,q) = ^ , , J+ — 

G-{p,po,q)^^—^ ^5±i;o±<?%+^ (85) 

irq{p + PQ + q) E+^q+p^Eq 



where G^{p,po,q) are the Laplace transforms of G^{x,t\xoO), corresponding to motion to the right 
y = x ~ Xq > and to the left y < 0, respectively (given that these will be later symmetrized). They 
precisely correspond to the exchange of p and po- We can now obtain the fully averaged diffusion front 
in our original model in the form: 

Gixo + y,t\xoO) = e{y)G+{xo + y,t\xoO) + e{^y)G-ixo + y,t\xoO) (86) 

as the sum of contributions of the G^" and (7^"+^ probabilities, respectively. The first one is obtained 
by setting po = — p, i.e. one has: 







dye-pyG+{xo + y,t\xoO) (87) 



2 



POO nL, pLi — Xq 

/ dLLe-P^L-^ dxo dye'PyG+{xo + y,t\xoO) 

Jq Jo J-xo 

^P^-P^P^-Tr EtE^ 

while the second one gives dyePyG~{xo + y,t\xoO) = p'^G~ {—p,p, p), i.e. the same function of p. 
Hence the final result is already symmetric. The Laplace inversion is thus immediate and one finds, 
using again (|3]): 

Gixo + y, t\xo, 0) = ( - _ - -)g(|;/|)e-^l^l (89) 

= P(xo,t|a;o,0)^G(|2/|/r2)e-l^l(''+*') (90) 
with r = Int, using our previous result (|5| for the average return probability. At strictly coinciding 



points the two Green's function must coincide, see Eq. (80 1 both being then proportional to the same 
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averaged trace over the same eigenenergies. Note however the divergence of G{X) ~ l/(-y7rX^/^) at 
small X hence g{\y\) ^ l/(^/7r|j/p/^) (probability of degeneracies being related to the return probability 
of a Brownian motion). Hence the result (90 1 is meant to hold in the scaling region y ~ and to 
break down for y — 0{1) since the true function g{y) should equal g{y = 0) = 1. At the end the result 
(90 1 is very similar to the one of pTj and [8 in the case p = 0, up to an exponential factor e"'''^', i.e. 
the absorbers act simply as an increase in the magnitude of the drift. 



10. Conclusion 



In summary we have computed various observables for the Sinai model in presence of a small density 
p of perfect absorbers, using the strong disorder RSRG method. We have confirmed the result of TH 
for the power law decay in time of the average return probability, and for the density of states of the 
corresponding Fokker-Planck and Schrodinger operators. We found that in presence of a small drift 
the decay remains a power law with a larger exponent. We have also computed the diffusion front 
and described the crossover from smaller times where the diffusion scale grows as L{t) Tint as 
in the absence of the absorbers, to the very large time limit when the survivors are found to be at 
a distance of order S^(t) = Tlni/y^p of their starting point, a new survival length scale. Remarkably 
their asymptotic distribution is a step function symmetric around the starting point, and even more 
remarkably it remains so in presence of a drift. Beyond this surviving zone around the starting point 
we found that the probability decays - in a scaling sense - also with a power law but with a larger 
exponent continuously increasing with the rescaled distance x/^{t). A qualitatively similar behaviour 
was found in the meeting and survival probability of two particles. The statistics of eigenstates was 
also studied using the RSRG and the connection with the methods used by Texier was emphasized. It 
was found that keeping only the lowest eigenstate in each interval between scatterers leads to the same 
leading large time result. From there a saddle point argument was presented to clarify the mechanism 
for power law decay in this model, as resulting from the competition of large regions free of absorbers 
and their anomalously long decay time. Finally, we computed the average Green's function of the 
associated Schrodinger problem, which is found to exhibit the same spatial decay as in the absence of 
absorbers, up to a global e"'''^' factor. 

There remains a host of quantities to compute for this problem using RSRG, as was done in Ref. 
pT| . It would also be interesting to develop other methods, such as the one used in p6] |27l |28] to 
compute diffusion fronts and first passage probabilities. There also remains fundamental questions, and 
let us list a few. Since a Sinai walker equilibrated for a time t — e^^'^ in a valley visits a site at potential 
U above the bottom of the valley proportionally to te~^/'^ , which is large for J7 < F, it is reasonable 
to expect that imperfect scatterers will renormalize to perfect ones. However a detailed study of this 
effect would be interesting. Other questions are to understand how general is the mechanism for the 
restoration of the symmetry in the presence of a drift observed here for some observables, what are 
the detailed properties of the eigenfunctions and the a.c. transport properties, whether the quenched 
versus annealed (in the probabilist sense) transition discussed in [25] also occurs here, and whether this 
model and some of the phenomena arising here can be extended to higher dimension, or to quantum 
models. 



Note added: after the completion of this manuscript we were communicated unpublished notes by 
A. Comtet, C. Texier and Y. Tourigny who obtained, via different methods and for a special value of 
the drift in the biased case, a similar decay exponent for the return probability as obtained here. 

Acknowlegments: I am grateful to C. Texier for a careful reading of the manuscript and useful 
remarks. I thank C. Hagendorf and A. Rosso for interesting discussions. 
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Appendix A. Survival of two particles 

We may ask to compute the probability that two particles distant by yo > be both still alive at time 
t. For simplicity we will compute the contribution, which we denote S{yo,t), to this probability which 
comes from the event when the two particles start in the same segment between absorbers. Hence this 
is a lower bound to the total probability. However, since when they start in different regions we can 
expect a faster decay, it already contains interesting information. 

Again, as in section [7] both initial conditions must be in the bulk region in Fig. [TJeach away from 
the edge bonds, i.e the absorbing zone. Hence it reads: 

S{yo, t)^p^ [ LdLe-P'^ I ^ < xa)e{xo + yo<L- £2^+2))^ (A.l) 
Jo Jo ^ 

where (..)^ denotes the average with respect to the fixed L measure ([ij). We have simply multiplied 
the probability that the leftmost particle happens to be on a segment of length L, with the uniform 
probability dxo/L that it is at position xq within dxo, and expressed the constraint. Introducing the 
Laplace transform w.r.t. yo and performing the integration ^ ^''^^ c?j/oe~™" ^''^^ ^° dxo one 
obtains: 

S{p,t) - / dy^e'^'y°S{yo,t) = -(ci>(p) _ ^$'(0) - $(0)) (A.2) 
Jo P 

^ n 



E-{p)E^{p) 
{p+p)E^{p + p)E+{p + p) 



where F = Tint, and S{y = 0, i) = S{t) the single particle survival probability computed in Sectionjsj 
In the symmetric case (up to an immaterial constant): 

<^{p)^pt^nh\T^p)—jj^—— (A.4) 
smh (1 y'p + p) 

Upon Laplace inversion one finds: 

<^{y) = pt'^nh'{V^)^ 2 (l-2^^m2^)e-(f^'"'+^)'' (A.5) 



m— — 00 



2.1 x-^ T _„2r 



pt&nh\T^p)—— '^"e-" -~py (A.6) 



^ n— — 00 



We now have: 



Vo 



S{y^, t) ^ 2 y^;: - pyo + / dz dz'^z') (A.7) 

smh(2Fyp) cosh^(F^) Jo Jo 

^ ^^P - Tr^T^m'^{2pya + 3) - 27r^m'^7/o _i^El^2^p)y„ ,^ 

Through an analysis very similar to Section [7] one finds the large time behaviour: 

S{yo, t) « iTlnt^pil - p^)t-2^^ (A.9) 



for < yo < (Tint)/ y/p, and a larger exponent for j/o > Tint/ ^fp, i.e. a decay S'(yoi t) ^ t 
with yo = yo^/{T\nt). 
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Appendix B. free diffusion 

For free diffusion (with unit diffusion coefRcient) with absorbing boundaries on [0,L] there are two 
dual formula for the diffusion front: 

PL{xt\xoQ) = y {e- '""^'"'''' - g- <-+-°+/"^'' )g(o ^ ^ ^ ^)^(o <xo<L) (B.l) 

v47rt 

* n— — oo 

9 cx) 

= 4 E sin(^) sin(^)e-'™'*/^'^(0 <x< L)e{Q < xo < L) (B.2) 
Ij Lj Lj 

Considering now the problem of Poissonian absorbers, one averages over the distribution of interval 
sizes L between absorbers and obtains the diffusion front averaged over (uniformly distributed) initial 
positions as: 

P{y,t) = p'j^ dLLe-P'^e{\y\<L) -^Pi^{x<y + y,t\xoO) (B.3) 



oo „c 



^^^-pL-.^™^*/^=((l _ M)cos(^) + ^sin(^)) (B.4) 
I L L m-K L 

From this expression, one obtains the average return to the origin probability, and its large time 
behavior; 

P{0,t)=p' V r dLe-^^--'-'^/^' « ?!^!;!^!p4/3,i/6e-3(f)^'%^^^t^^^ (B.5) 

retaining m = 1 and using the saddle point at i = 2^/^7r^/^t^/^/9~^/^ and computing the fluctuations 
around it. One also obtains the survival probability, and its large time behaviour as: 

S{t) = r dyP{y, t) = p2 r dLLe-P"^ f ^ ^ dxP^x, t\xo, 0) (B.6) 

J-oo Jo Jo ^ Jo 

= p'y rfiL^(l-(-ine-''^-^-^*/^^«4^piV2e-3(f)=^W/3 

m=l " * 



using again the saddle point and m = 1. At large time the full diffusion front is also dominated by 
m = 1: 

P(y,i)«yp2 / rf^e-p.L-.Wi-^)((i_ ^eos(^) + -sin(3)) (B.8) 
Ji L L 1^ L 

with L = L/y. Upon a similar saddle point analysis one finds that the diffusion front expressed in the 
variable: 

y = y/Ut) , Co(t) = 2^3^2/3^-1/3^1/3 (B.9) 

takes a simple scaling form: 

P{y^t)dy = 5(i)^((l - |y|)cos(7r|y|) + ^ sin(7r|y|))0(l - \y\)dy (B.IO) 

Hence the survival zone has also sharp edges at y = ±1, though this time it is continuous at the 
edge (the rescaled probability vanishes as ^ (f — lyl)"* at the edge). By the same mechanism as 
described in the text, the decay of the average front outside this region is as exp{—C{y)t^^^), with 
C{y) = |C(l)(y + p^), i.e. the amplitude in the stretched exponential grows continuously with y > 1. 



Sinai model in presence of dilute absorbers 



21 



References 

[1] Ya. G. Sinai, Theory of Prob. and Appl. 27(2), 247 (1982). 

[2] J. -P. Bouchaud, A. Comtet, A. Georges, and P. Le Doussal, Europhys. Lett. 3 653 (1987) and Ann. Pliys. (N.Y.) 
201, 285341 (1990). 

[3] David K. Lubcnsky, David R. Nelson, arXiv:cond-mat/00 04423l Phys. Rev. Lett. 85, 1572 (2000) and |arXiv:cond-| 

mat/0107423 D. R. Nelson, arXiv:cond-mat/0309559 
[4] I. Aranson, L.Tsimring, V. Vinokur, adap-org/9702002 

[5] D. Dufresne, Scand. Act. J. (1990) 39, H. Geman and M. Yor, Math. Fin. 3 (1993) 349. 

[6] L M. Lifshits, S. A. Gredeskul, and L. A. Pastur, Introduction to the theory of disordered systems, John Wiley 
and Sons, 1988, A. A. Gogolin and V. I. Melnikov, Sov. Phys. JETP 46, 369 (1977). A. A. GogoUn, Phys. Rep. 
86(1), 153 (1982) 

[7] C. Texier, J. Phys. A: Math . Gen. 33, 6095 (2000^^ 

[8] L. Balents, M. P. A. Fisher, arXiv:cond-mat /9706069I 

[9] H. Kesten, M. Koslov and F. Spitzer, Compositio Math. 30 (1975) 145, H. Kesten, Physica 138 A (1986) 299. 
[10] A.O. Golosov, Comm. Math. Phys. 92 (1984) 491 

[11] P. Le Doussal, C. Monthus, and D. S. Fisher, condmat/9811300, Phys. Rev. E 59(5), 4795 (1999). D. Fisher, P. 

Le Doussal and C. Monthus, Phys. Rev. Lett. 80 (1998) 3539. 
[12] D. S. Fisher Phys. Rev. B 50, 3799 (1994) and D. S. Fisher Phy s. Rev. B 51, 6411-6461 (1995). 
[13] For a general review see: F. Igloi, C. Monthus, a rXivxond-mat /0502448| Physics Reports 412, 277-431, (2005) 
[14] G. Texier and C. Hagendorf arXiv:0902.2698 

[15] G. Hagendorf and G. Texier, J. Phys. A: Math. Theor. 41, 405302 (2008). 
[16] D.S. Fisher and A. P. Young, Phys. Rev. B 58, 9131 (1998) 
[17] G. Monthus, arXiv:cond-mat/0309029 , Phys. Rev. B 69, 054431 (2004). 
[18] G. Monthus, arXiv:cond-mat/0212212 , Phys. Rev. E 67, 046109 (2003). 

[19] see Sec. VIILA in G. Monthus and P. Le Doussal, arXiv:cond-mat/0202295 , Phys. Rev. E 65 (2002) 66129. 

[20] A. Gomtet, G. Monthus, M. Yor, arXiv:cond-mat/9601014 J. Appl. Proba. 35 (1998) 255. 

[21] see e.g. Ref. [2l 1111 for discussion of this modcl-dcpcndent crossover at very short times, not studied here. 

[22] at larger p, it remains to be understood if the coincidence of prefactor with TH is accidental or not. 

[23] this amounts to neglect relaxation on times scales tj = 1/Ej w.r.t. ti with j > i. This is the standard 
argument of RSRG being asymptotically exact: this neglect can only produce an uncertainty in of order 
rin(ti ±tj) = Tlnti -|-ln(l ±e~('""~^j'/^) = F^ -|- 0(e~(^'~^j)/^). For decimations which occur in the same 
region of space (i.e. overlapping eigenfunctions) Fi — Tj is of the same order as F; hence the error is exponentially 
small in F at large F. 

[24] D. S. Dean, G. Sire, J. Sopik, 'arXiv:cond-m at/0604456t Phys. Rev. E 73, 066130 (2006) 

[25] P. Le Doussal and G. Schehr, in preparation. 

[26] D. Dean and A. Gomtet, J. Phys. A31 (1998) 8595 

[27] S. Majumdar and A. Gomtet, Phys. Rev. E 66 (2002) 061105. 

[28] G. Sire, arXiv:cond-mat/9902223 Phys. Rev. E 60, 1464 (1999). 

[29] G.Ben Arous, S. Molchanov and A. F. Ramirez, arXiv:math/0501107v 3, The Annals of Probability 2005, Vol. 33, 
No. 6, 2149. 



